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1 Introduction 



Abstract We study entanglement dynamics in the presence of correlated environmental 
noises. Specifically, we investigate the quantum entanglement dynamics of two spins in the 
' (— I ' , presence of correlated classical white noises, deriving Markov master equation and obtain- 

^ ^' ing explicit solutions for several interesting classes of initial states including Bell states and 

X form density matrices. We show how entanglement can be enhanced or reduced by the 
P5 ■ correlation between the two participating noises. 

■ 

^ ' Keywords Entanglement dynamics • Correlated noises • Two-qubit model • Entanglement 

sudden death 

PACS 03.67.-a ■ 03.65.Yz ■ 03.65.Ud ■ 03.67.Mn. 
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' Many proposed applications in quantum computing (T|, quantum communication (2||, and 

, quantum cryptography |3| revolve around harnessing the inherent correlation between quan- 

' tum particles, called entanglement |4]. Although quantum mechanics dictates that these co- 

^\ \ herence effects are intrinsic in certain systems, even when the atoms or particles are non- 

■ local, there is an overall weakening due to coupling to noisy environments f5l that eventually 

leads to the fast decay of entanglement [6j.7j^,9,,_10J, as in the case of amplitude or phase 
noise, or even the sudden death of entanglement (ESD) in the worst scenarios llllll2lfT3]ll4l 
' [T5l. As such, the study of the controlled entanglement dynamics is of much importance to 

the prospects of maintaining quantum information 1 16, 17, 18, 19, 20, 21, 22, 23 1. Moreover, 
model systems that theoretically exhibit the rebirth of entanglement have been proposed and 
discussed in several cases I24l|25|[26ll27| . 

In this paper, we consider a model system consisting of two uncoupled qubits A and B 
interacting with stochastic fields, /^(f) and fsit), respectively, shown in Fig.l. The setup 
of the system is quite similar to the two-qubit local dephasing channel of 1 8 1 in which the 
qubits were found to disentangle in a shorter time than their individual atom-field local de- 
phasing times. On the other hand, if the two qubits are coupled to a common environment. 
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then entanglement was shown to be preserved for a class of initial states living in a subspace 
called decoherence-free subspace (DFS) |28,29]. By introducing a correlation between the 
two noisy sources, we can show that the entanglement between the two qubits can be en- 
hanced or reduced by properly choosing the correlation of the two participating external 
noises. 

The paper is organized as follows: We present the specifics of the model system in 
Section II, which leads to deriving the master equation in the Markov regimes. The entan- 
glement evolution under the correlated noisy sources are studied in Section III. We have 
shown that dependent on initial states the correlation between two external noises can either 
enhance or reduce the existing entanglement of the qubit systems. We conclude in Section 
IV. 



2 Qubit Model and the Master Equation 

Our model consists of two identical, separated qubits (two-level atoms, spins, excitons, etc.) 
each having transition frequency w and each coupled to a separate classical random field, 
/^(?) and /^(f) respectively. The total Hamiltonian for this system is (setting h= \) 

//tot = '^{at + a^)+fA{t)o^ + fB{t)a^. (D 

This Hamiltonian is familiar in condensed matter theory as an extension of the spin-boson 
model in the semiclassical regime [301. The classical stochastic fields are assumed to obey 
the following correlation relations in the Markov approximation; 

M[fA{t)]=M[fB{l)]=0, (2) 
M[fA(t)fA{s)] = YA8{t-s), (3) 

M[fB(t)fB{s)] = YB8{t-s), (4) 
M[fA{t)fB{s)]=rd{t-s), (5) 

where M[-] denotes the ensemble average over the classical stochastic fields and for simplic- 
ity we assume Ya = Yb = Y- Obviously, F determines the respective correlation properties. 










fA(t) 




fB(t) 









Fig. 1 The model is composed of two uncoupled qubits A and B individually coupled to respective stochastic 
fields, /A(f) and /s(f), which are then correlated. 
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The effects of the correlation between the two fields /4(r) and /^(r), M[fA{t)fB{s)], will be 
revealed in the study of entanglement between the two qubits. 

There has been a lot of work dedicated to quantum entanglement dynamics in the pres- 
ence of an environmental noise, either in classical or quantum regimes (For some recent 
progress, e.g., see I31ll32ll33]|34ll35ll3"6ll37ll38|[39j ). Clearly, having two separated qubits in- 
dividually coupled to respective external fields would reveal the same Lindblad dynamics 
for each qubit. That is, individual couplings always cause the irreversible decay of entan- 
glement. However, when there is a correlation between the two classical stochastic fields, 
we expect to find new entanglement effects between the two qubits. It can be shown that the 
master equation governing the dynamics of two qubits under the influence of two correlated 
noises can be derived from the corresponding stochastic Schrcidinger equation 1401 , 

p = -i[Hs,p] (6) 
-27a (P - o^pa^) - 2yb{p - o^pa^) 
-4r{a^a^p + pa^a^ - a^pa^ - a^pa^). 

where Hs = ^a{ai^ + a?) is the Hamiltonian of the two qubits. Upon expanding and 
af in terms of CT^ *, we may obtain a more solvable form to the master equation: 

p = -i[Hs,p] -YAia^a^p-a^pa^-a^pa^ + pa^a^) (7) 
-Ysia'^a'^p - a'^pal - a'^pa'^ + pa'^a'^) 
-rio^o'^p - a'^pal - o^pa'^ + pa'^a^) 
-r{a1a^p - a^pal - alpa^ + pa^al)+H.C. 

At a first glance, Eq. ^ or Eq. ^ immediately displays the dynamics of qubit A in 
correlation with field fA{t), dynamics of qubit B while interacting with field and the 
cross terms due to the correlation between the two fields. Without the involvement of F, 
the master equation would be the sum of two familiar Lindblad master equations for two 
qubits, respectively. In that case, the two qubits would evolve separately throughout time and 
entanglement between the two qubits will deteriorate with time. The correlation between two 
classical stochastic fields characterized by F is expected to affect entanglement evolution in 
two different ways dependent on the initial states. We will show for some initial states that 
entanglement can be significantly enhanced by adjusting the cross correlation between the 
two fields. However, for some other initial states, the cross-correlation works like a catalyst, 
which may accelerate the decay of entanglement. 



3 Modulated Entanglement Evolution 

Modulated entanglement evolution will of course depend on the state the system was origi- 
nally in. If the main goal is to try to maintain or improve entanglement, a good initial state 
would be that which has maximum entanglement, the Bell State, as is presented in cases 
13.11 and 13.21 It is also interesting to see how a general X type matrix behaves under these 
circumstances, as in case 13.31 It is important to note that y and F cannot be arbitrary and 
must be chosen in a way that preserves positivity of the density matrix (F < y). 
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3.1 Bell State (| ++), 



-» 



One of the Bell States is described by the pure state vector 



1 



W = ^{K+> + |- 

which has the following density operator representation: 






VI/ 



^(1001)4 



-)}, 



/I 00 1\ 
0000 
0000 

VI 00 1/ 



(8) 



(9) 



Inserting p as the initial reduced density matrix, p(0), fuels the following master equa- 
tion solutions: 



Pn[t) 
P22(r) 



1 

8k- 



„-(6y+rc)/r 



e ■ [27(1 - e^"") + K{2e^^^+'''^' + e^"" + 1)], 



8k: ^ 



+ k:(2, 



,(67+)c)r _ IKI _ 



1)], 



Pi4(r) 



where we have defined k: = 1/47^ + 32F^. It should be noted that this special initial condi- 
tion leads to the simple relations: 

P33(f) =P22(r) P44(f)=Pll(0 P32(0=P23(0 P4l{t) = Pu{t) , 

Pl2(0 = Pl3(0 = P24(0 = P34(0 = 0, 
P2l(0 = P3l(0 =P42(0 =P43(0 =0. 

Because the off-diagnal elements initially at zero will remain at zero for all time, the 
time-dependent reduced density operator will be of the X form with only 4 independent, 
real terms: 

fpu{t) pu{i)\ 

P22(f) P23(0 
P23(0 P22(0 

VPi4(0 pu{t)J 

With a solution to the master equation in hand, it is easy to calculate the Concurrence 
1411 . a measurement of entanglement ranging between and 1 evaluated as 



Pit) 



(10) 



C(p) = 2max{0, V^i " V^2 - Vh - V^4}, 
where A; represent the eigenvalues of the matrix 

p = p(a^®a,^)p*{a^®a^), 
in descending order. For this p(f), the concurrence throughout time will be 



(11) 



C(p(0) = 2max{0, IpuWI - VP22WP33(0}> 



(12) 



Fig. 2 Concurrence calculated for various values of F. At times later in the decay of the concurrence, entan- 
glement can be built up by increasing the correlation between the fields. 



and is plotted against jt in Fig. 2 for various values of T . 

For all cases of y, when F = 0, the system exhibits the sudden death of entanglement 
well-known to the local-dephasing channel 1 8 1 . However, by turning on the correlation F 
between the stochastic fields, the concurrence curve moves vertically upward, denoting an 
increase in the measurement of entanglement. Clearly, we displayed controlled entangle- 
ment evolution via correlated environmental noises. We also are faced with the fact that 
increasing the correlation y of each atom to their respective fields will also induce a faster 
sudden death of entanglement. By turning on the correlation between the fields to its maxi- 
mum value, F = y, entanglement is maximally enhanced, as shown in Fig. 2. 



3.2 Bell State(|+-),|--F)) 

Let's now solve the master equation using the other form of the maximally coherent Bell 
State: 



lc&o) = ^{| + -) + !-+)}. 



(13) 



The initial reduced density matrix for this state vector is: 




/0000\ 
110 
110 

Voooo/ 



(14) 



which allows us to solve for the time dependent density matrix: 



Pl2(f) = Pl3(0 = P24(0 = P34(0 = Pl4(0 = 0, 
P2l(0 = P3l(0 =P42(0 =P43(f) =P4l(0 =0. 



(15) 



Fig. 3 Concurrence calculated for various values of F with 7=1. For this initial state, entanglement can 
only be reduced by turning on the correlation between the fields. 



Using the same definitions as before 



Pii(0 

P22(0 
P23(0 



1 

iff 
1 

1 

4k 
1 

4^' 



^-(6y+K), [_ (2y + 8r) ( 1 - e^-^' ) - 1 + e^"' " 2e'^^''+''^' )] , 
-(67+^)rj(2y+8r)(l - e^"') + k{1 + e^"' +26^^^+"^')], 



-{6Y+K)t 



[- (4y - 8r ) ( 1 - ) + 2 k:( 1 + e"" ) 



P33(0=P22(0 P44(0=Pll(0 P32(0=P23(0- 

Making tfie density operator of a particular X form with 3 independent, real variables 
throughout all time t: 



Pit) 



/pii(r) \ 

P22(0 P23(0 
P23(r) P22(0 

V pii(07 



(16) 



The concurrence of this density matrix l |16l l is given by: 



C(p(f)) = 2max{0, IP23WI - VPuJtJp^}, 



(17) 



as plotted in Fig. 3. 

This figure displays quite the opposite of that in the previous section. Introducing a 
correlation between the noise fields, F, actually lowers the concurrence curve and causes a 
faster decay of entanglement. This is true for all values of y, displaying the impossibility for 
entanglement enhancement for this initial state evolution. This brings forth a very interesting 
characteristic of entanglement, that even in the Markovian regime where there is no memory 
of previous times in the system, it is still very sensitive to which initial state is being used. 
Even though all Bell States produce maximum entanglement, one has the prospect for the 
rebirth of entanglement and the other does not. 



7 



3.3 X state 



We can now use a more general approach by utilizing an initial state in the X form which 
includes the Bell states and Werner states as special cases |42 |: 



/ fl w* \ 

bz* 

z c ' 
VwO d J 

where a + b + c + d=l.ln this case, the concurrence is 

C{Pa) = 2max{0, |z| — \fad, \w\ — \fhc\, 

which allows us to account for a large range of initial entanglement conditions. 
The solution to the master equation is: 



(18) 



(19) 



■ + 



■ + 



e-6y/ 
Ak 


[(«- 


h- 


c + c;)K:coshK:f + (16rz- 


-ly{a- 


h- 


c + d)) sinhiff 


4k- 


[(«- 


h- 


c + (f cosh (Cf + ( 1 6rz - 


2r(«- 


h- 


c + J))sinhK'r' 


Ak 


[(«- 


h- 


c + (i) (f cosh (Cf + ( 1 6rz - 


2r(«- 


b- 


c + (i))sinh(c;' 


e-6r' 
Ak 


[(«- 


h- 


c + (i)K-coshK-f + (16rz- 


-2y(a- 


h- 


c^d)) sinhiff 



PiiW = ^ — — 

Piiit) = ^ 

P33{t) = ^ 

2(d-a)e-'^y 
P44(0 = ^ '— 

Pai') = e"''^'[zcoshK-f + -((a-Z7-c + J)r + 7z)sinhK:r] =p32(0i 

Pl4(0 = we-'^^' = P4i{t). 

The density operator describing the mixed states of the system will now remain in the X 
form throughout time with 6 independent, real terms 



Pit) 



/pu{t) pi4(r)\ 

P22(f) P23(0 
P23(0 P33(0 
\ Pl4 P44(0y 



(20) 



such that the concurrence is calculated as 

C(p(f)) = 2max{0, |p23(0l - VPn (0P44(r), IPmWI - \/p22(0P33(f)}- (21) 

The X form initial density matrix displays a more general range of control over the 
modulation of entanglement. As shown below, the correlation may cause entanglement to 
decay faster initially, but it can compensate the loss of entanglement later. To see this, we 
choose specific values of the initial X form density matrix 



P(0)=3 



110 
110 

Viooi/ 



(22) 



and plot the concurrence through time in Fig. 4. 
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0.35 




Fig. 4 Concun'ence plotted over yt for various values of F. At first entanglement is reduced due to correlated 
noise but later begins to enhance it. 



The results demonstrate that the entanglement can be significantly enhanced by the 
cross-correlation at late times even though the correlation was detrimental to entanglement 
at early times. More explicitly, it can be shown that at early time |p23 (0 1 ~ i/pii (r)p44(f) is 
the dominant term in the calculation of concurrence, resulting in the degradation of entan- 
glement due to the cross-correlation. However, at a certain point \p\i,(t) \ — \J Pn.[t)pyi(t) 
begins to dominate, causing the correlation to enhance the entanglement. This demonstrates 
our ability to improve entanglement for a wide range of initial states and time scales. 



4 Conclusion 

In summary, we studied a system of two separated qubits each coupled to a stochastic field 
which can then be correlated, opening the option for the enhancement or reduction of entan- 
glement between the qubits. By solving the master equation for the qubit dynamics under 
various initial states and viewing the concurrence as a function of time, the effects of the 
correlation between the fields, F, became imminent. Of our three cases for initial states, the 
Bell State, Eq. showed a potential for the regeneration of entanglement. The amount 
of entanglement this particular case can achieve is dependent on the degree of correlation 
between the two noises. In the second case we considered the Bell State of Eq. i ll4b for 
which correlation F did not cause the entanglement between the atoms to be regenerated. In 
fact, it instead reduced the entanglement and caused it to decay at an even faster rate. This 
illuminates a selectivity that entanglement has toward the initial state of the system, even 
under the memoryless Markov approximation. It is important then to look at the general X 
form density matrix of Section [331 allowing us to account for various initial entanglement 
conditions of the system. This case presented the ability to modulate entanglement in either 
way. The range of maximum enhancement was variable depending on the parameters of the 
initial matrix, giving us a wide range of control over the improvement of entanglement. 

In general, many characteristic properties of entanglement are still unknown. In this 
model system, we have created a bridge from one qubit to the other through the correlation 
of the stochastic fields, allowing for entanglement to be modulated. This correlated noise 
approach can easily be applied to multipartite systems, which in turn might provide an even 
higher enhancement of entanglement |43 ,44 ,45 ,46 1. Finally, it may be worth noting that the 
relationship between the classical noise model presented here and the fully quantized models 
discussed in I17III8I is an interesting problem that will be addressed in future publications. 
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